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ABSTRACT 

We have constructed MOCASSIN photoionization plus dust radiative transfer models for the Crab Nebula core¬ 
collapse supernova (CCSN) remnant, using either smooth or clumped mass distributions, in order to determine 
the chemical composition and masses of the nebular gas and dust. We computed models for several different 
geometries suggested for the nebular matter distribution but found that the observed gas and dust spectra are 
relatively insensitive to these geometries, being determined mainly by the spectrum of the pulsar wind nebula 
which ionizes and heats the nebula. Smooth distribution models are ruled out since they require 16-49 M 0 of 
gas to fit the integrated optical nebular line fluxes, whereas our clumped models require 7.0 M 0 of gas. A global 
gas-phase C/O ratio of 1.65 by number is derived, along with a He/H number ratio of 1.85, neither of which can 
be matched by current CCSN yield predictions. A carbonaceous dust composition is favoured by the observed 
gas-phase C/O ratio: amorphous carbon clumped model fits to the Crab’s Herschel and Spitzer infrared spectral 
energy distribution imply the presence of 0.18-0.27 M 0 of dust, corresponding to a gas to dust mass ratio of 26- 
39. Mixed dust chemistry models can also be accommodated, comprising 0.11-0.13 M 0 of amorphous carbon 
and 0.39-0.47 M 0 of silicates. Power-law grain size distributions with mass distributions that are weighted 
towards the largest grain radii are derived, favouring their longer-term survival when they eventually interact 
with the interstellar medium. The total mass of gas plus dust in the Crab Nebula is 7.2±0.5 M 0 , consistent 
with a progenitor star mass of ~ 9 M 0 . 

Subject headings: ISM: supernova remnants: individual (Crab Nebula); circumstellar matter 


1. INTRODUCTION 

Evolved stars, in particular AGB stars, have long been con¬ 
sidered as significant contributors to the dust found in the in¬ 
terstellar media (ISMs) of galaxies. However, recent quan¬ 
titative determinations of AGB star dust injection rates into 
the ISMs of nearby galaxies such as the LMC have found sig¬ 
nificant shortfalls compared to current estimates for the re¬ 
quired replenishment rates, e.g. Matsuura et al. ( |2009| l; Boyer 
et al.| < 2011 1 ; Matsuura et al. (201 3|>. Influenced in particular 


by discoveries of very large dust masses in some high red- 
shift galaxies emitting less than a billion years after the Big 
Bang ( Omont et al.||200T[ Carillr et al.||200l[ Bertoldi et al. 
|2003| l, the potential contribution of core-collapse supernovae 
(CCSNe) to ISM dust budgets has also been investigated in¬ 
tensively in recent years. To significantly influence the dust 
budgets of galaxies, ejecta dust masses of at least 0.1 M 0 
per supernova have been judged necessary (|Morgan & Ed 


munds|2003[|Dwek et al.|2007[|Michalowski et al. |20 1 01 |Gall 

et al.||2011) . While some CCSN dust formation modellers 


have predicted that such masses of dust should form, others 


have not (Kozasa et al. 1991 ;|Todini & Ferrara 2001 

Nozawa 

et al. 

2007 

Bianchi & Schneider 2007; 

Sarangi & C 

aerchneff 

2013 

2014 

). Observational determinations of how much dust 


has formed in the ejecta of CCSNe are therefore key. 


Starting with SN 1987A (e.g. 

Wooden et al.| 1993; Bouchet 

& Danziger 1993; Ercolano et a 

. 2007a) and continuing with 

Spitzer stud 

ies of CCSNe (e.g. Sugerman et al. 

2006 ;JMeikle 

et al. 

2007 

Andrews et al. 2011. Kotak et al. 

2009l]Fabbri 

et al. 

2011 

Meikle et al. 2011), mid-infrared observations 


during the first 3-4 years after outburst typically measured 
no more than ~ 1 0 -3 M 0 of newly formed 200-450 K warm 
dust in the ejecta, well short of the quantities required for 
CCSNe to significantly influence galaxy dust budgets. Re- 
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cently however the Herschel Space Observatory has detected 
large masses of much cooler dust within several young core¬ 
collapse supernova remnants (SNRs). Barlow et al. ( 2010} 
measured 0.075 M 0 of cool ~35 K dust emitting at wave¬ 
lengths longwards of 70 pm in the Cassiopeia A SNR, which 
together with the 0.025 M 0 of warm dust measured by Spitzer 
to be emitting shortwards of 70 pm ( jRho et al.|2008 1 implied 
a total of 0.10 M 0 of new d ust within this 340-yea r old SNR. 
Following the discovery by Matsuura et al. ( 2011 [ > with Her¬ 
schel of 0.4-0.7 M 0 of cold dust in the then 23-year old rem¬ 
nant of SN 1987A, high angular resolution ALMA observa¬ 
tions at 440 and 870 pm ( |Indebetouw et al.[2014[ ) confirmed 
that the cold dust was located in the ejecta, with a mass of 0.5- 
0.8 M 0 ([Matsuura et a l.|2014| . This implied a large increase 
in the ejecta dust mass during the more than 20 years that had 
elapsed since the mid-IR observations that had detected less 
than ~10 -3 M 0 of warm dust (Wesson et al. 2015). From 
Herschel ob servations of the 960-year old Crab Nebula SNR, 
Gomez et al. (j2012| deduced the presence of 0.12 M 0 of 
amorphous carbon or 0.24 M 0 of silicates, much larger than 
the ~ 3 x 10 1 M 0 of warm dust that had been derived from 
shorte r wavelengt h Spitzer observa t ions (|Temrm et al.||2006[ 
2012|l. However, Temim & Dwek (201 3]> subsequently pre¬ 


sented radiative transfer modelling of the Spitzer and Herschel 
observations of the Crab Nebula, assumimg a central point 
heating source, to obtain lower dust mass estimates, namely 
0.02-0.04 M 0 of amorphous carbon, or 0.13 M 0 of silicates. 

The Herschel observations of the Cas A, SN 1987A and 
Crab supernova remnants that have been summarised above 
have refocused attention on the potentially significant contri¬ 
bution that core-collapse supernovae can make to interstellar 
dust budgets. Given the importance of an accurate dust mass 
estimate for the Crab Nebula, we have constructed a number 
of gas+dust radiative transfer models for the nebula that use 
the diffuse radiation field of the pulsar wind nebula (PWN) 
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Fig. 1.— Composite image of the Crab Nebula, obtained by combining a Hubble Space Telescope optical emission line image (blue-white) with a Herschel 
Space Observatory 70-/im dust emission image (red), showing the emitting dust to be closely aligned with the optical knots and filaments. The image is 
5.65 arcmin on a side; north is up and east is to the left. Credits: Oli Usher (UCL); Herschel Space Observatory, Hubble Space Telescope: ESA, NASA. 


along with realistic nebular geometries and density distribu¬ 
tions. In these models, dust grains with a range of composi¬ 
tions and size distributions are immersed in nebular gas out¬ 
side the PWN, with the gas either (a) smoothly distributed, 
or (b) within clumps that mimic the Crab’s highly filamentary 
structure. We present first our fits to the integrated optical 
emission line fluxes measured by |Smith[ ( |2003| ) for the Crab 
Nebula, yielding gas-phase elemental abundances and masses 
in the nebula. We then present the results from our modelling 
of the infrared spectral energy distribution using several dif¬ 
ferent potential grain species, and compare our derived dust 
masses with previously published dust mass estimates. 

2. INPUT PARAMETERS FOR THE GAS AND DUST MODELS 


The Crab Nebula is the remnant of a supernova that was 


recorded in 1054. A distance of 2 kpc is often adopted (Trim¬ 
ble: 1968 I. It is one of the best-studied objects in the sky, hav¬ 


ing been observed at all wavelengths from 7 -rays to the radio. 
It has been suggested to have resulted from a Type Iln-P core¬ 
collapse explosion of a progenitor star whose initial mass was 
~10 M 0 (|Smith||2013|>. The nebula is rare amongst SNRs 
in not being collisionally ionized but is instead photoionized 
by synchrotron radiation from the pulsar wind nebula at the 

center of the remnant ( |Hester|2008|). _ 

We have used mocassin ( |Ercolano et aT7| |20Q3[ |2005[ 


2008), a 3D photoionization and dust radiative transfer code 


which allows for arbitrary gas and dust geometries and den¬ 
sity distributions, diffuse radiation fields and multiple point 
input radiation sources with user-specified spectra, and mul¬ 
tiple dust grain species having user-specified grain size dis¬ 
tributions. MOCASSIN self-consistently solves the equations 
of radiative transfer to determine within each cell the degree 
of ionisation and the gas and dust temperatures, along with 
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TABLE 1 

Model parameters 



Model I 

Model ll 

Model III 

Model IV 

Model V 

Model vi 

Density distribution 

Smooth 

Smooth 

Smooth 

Clumped 

Clumped 

Clumped 

Total dimensions 

4.0x2.9 pc 

4.0x2.9 pc 

4.0x2.9 pc 

4.0x2.9 pc 

4.0X2.9 pc 

4.Ox2.9 pc 

Inner axes 

1.lx 1.1 pc 

2.1 x 1.4 pc 

2.3 x 1.7 pc 

3.Ox 2.0 pc 

2.3 x 1.7 pc 

2.3 x 1.7 pc 

H-density 

1400 cm -3 

775 cm" 3 

675 cm -3 

1700 cm' 3 

1900 cm -3 

1900 cm -3 

Radius of each clump 




0.037 pc 

0.037 pc 

0.037 pc 

Final gas-phase abundances, by number 

Hydrogen 

1.00 

1.0 

1.00 

1.00 

1.00 

1.00 

Helium 

1.83 

1.90 

1.90 

1.83 

1.85 

1.85 

Carbon 

9.7xl0" 3 

9.3x 10 -3 

9.3xl0 -3 

9.7xl0 -3 

1.02x 10 -2 

1.02x 10 -2 

Nitrogen 

2.5x 10 -4 

1.5x 10 -4 

1.5xl0 -4 

2.5x 10 -4 

2.5xl0 -4 

2.5x 10 -4 

Oxygen 

7.2x 10 -3 

8.0x 10 -3 

7.0x 10 -3 

7.2x 10 -3 

6.2xl0 -3 

6.2xl0 -3 

Neon 

2.0x 10 -3 

4.5x 10 -4 

4.5xl0 -4 

2.0x 10 -3 

4.9xl0 -3 

4.9x 10 -3 

Sulphur 

4.0x 10 -4 

5.0x 10 -5 

5.0x 10 -5 

4.0x 10 -4 

4.0xl0 -5 

4.0x 10 -5 

Argon 

5.0x 10 -5 

6.0xl0- 6 

5.0x 10 -5 

5.0x 10 -5 

l.OxlO -5 

l.Ox 10 -5 


the overall line and continuum output spectrum from X-ray 
to sub millimeter (submm) wavelengths of the region being 
modelled. We used MOCASSIN 2 .02.70 to fit the Crab’s ob- 
served infrared and submm SEP (fTemim et al.|2006[ |Gomez| 


et al.|201~2] Planck Collaboration|201 1 1 , along with the inte¬ 

grated optical nebular emission line fluxes measured by [Smith] 
(2003). 


was 

-l 


2.1. The input radiation field and the nebular geometry 

The adopted overall geometry for the nebula was an ellip¬ 
soid with a major axis diameter of 4.0 pc and a minor axis di¬ 
ameter of 2.9 pc (Hester 2008| l. The synchrotron-emitting pul¬ 
sar wind nebula permeates this volume, which is also partly 
occupied by the gas corresponding to the observed clumps 
and filaments. 

The pulsar wind nebula’s synchrotron spectrum from 
0.36 nm to 1 m that was used for the modelling was a digitized 
version of the spectrum plotted by Hester ( 2008| >. The level 
of the submillimeter part of the input spectrum needed to be 
lowered slightly in order to be consistent with recent Planck 
observations (Planck Collaboratio n[2011) . The spectrum 
scaled to have an integrated luminosity of 1.3 xlO 38 erg s 
(Hester p008| . The angular extent of the synchrotron emis¬ 
sion from the pulsar wind nebula appears to be a function of 
frequency, with the radio emission extending throughout the 
4.0x2.9 pc ellipsoidal nebula (Hester |2008| >, while at X-ray 
wavelengths the PWN has a diameter of ~ 1 pc (Hester et ak 
[ 2002 ) . 

To investigate the effects of different distributions of gas 
and dust within the nebula, several shell and PWN geometries 
were therefore investigated. Table [I] summarises some of the 
parameters used for the nebular models described below. 

I. a smooth shell distribution, with the gas and dust located 
at a radius of 0.55 pc in a 0.1 pc thick shell (i.e. both 
inner axes 1.1 pc in length), with the PWN diffuse field 
radiation field emitting uniformly from within the inner 
nebular radius of LT pc. This shell geometry was ar¬ 
gued for by Cadez et al. (2004) based on their multi- 
slit s pectroscopy and was adopted by |Temim & Dwek| 
( 2013| > for their dust modelling. A shell hydrogen den¬ 
sity of 1400 cm -3 was found to be needed to match the 
total (dereddened) H/3 flux fom the nebula. 

II. a smooth distribution of gas and dust in a shell with inner 


axis diameters of 2.lx 1.4 pc that extends to the outer 
nebular boundaries, immersed in a diffuse PWN radi¬ 
ation source that also extends to the outer boundaries. 


This corr esponds to the geometry discussed by Davidson 
& Fesen (1985). A shell hydrogen density of 775 cm _J 
matched the total nebular H/3 flux. 

III. a smooth gas and dust distribution in a shell that has inner 
axes of 2.3x 1.7 pc (Lawrence et a l.| 1995 ), extending all 
the way to the outer 4.0x2.9 pc limits of the nebula, as 
does the PWN diffuse radiation source. A shell hydro¬ 
gen density of 675 cm -3 was found to match the total 
nebular H/3 flux. 

IV. a clumped shell distribution that has inner axis diame¬ 
ters of 3.0x2.0 pc and extending to the 4.0x2.9 pc outer 
nebular edges, but in this case with the diffuse radiation 
source located entirely inside the inner axes of the shell. 
The degree of clumping is determined by fitting the op¬ 
tical line strengths. A clump filling factor of 0.10 and a 
clump H-density of 1700 cm -3 were found to be needed. 
The clumps are 0.037 pc (3.8 arcsec) in radius. The num¬ 
ber of clumps decreases with nebular radius as r -2 . The 
clumps are modelled using sub-grids, as described by Er- 
colano et al.| ( |2007b) . 


V. a clumped shell distribution where the gas and dust 
cl umps s tart at inner axis diameters of 2.3 x 1.7 pc 


(Lawrence et al. 1995), and with an r 2 distribution of 


clumps that extends to the 4.0x2.9 pc outer boundaries 
of the nebula, with a volume filling factor of 0.10. The 
PWN radiation field is a diffuse source emitting uni¬ 
formly within a 1 .lx 1.1 pc diameter sphere at the centre 
of the nebula. For 0.037-pc radius clumps, a H-density 
of 1900 cm -3 within the clumps was found to match the 
total H/3 flux from the nebula. 

VI. a clumped shell distribution where the gas and dust 
cl umps sta r t at inner axis diameters of 2.3 x 1.7 pc 


(Lawrence et al. 1995), and with an r 2 distribution of 


clumps that extends to the 4.0 x 2.9 pc outer boundaries 
of the nebula, with a volume filling factor of 0.10. The 
clumps are immersed in the PWN radiation field emit¬ 
ted from the entire volume of the nebula. For 0.037- 
pc radius clumps, a H-density of 1900 cm -3 within the 
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clumps was found to match the total H/3 flux from the 
nebula. 


Our preferred geometries are clumped Models V and VI. A 
clumped version of smooth Model I could not be constructed: 
its shell is only 0.1 pc thick and already required a relatively 
high H-density of 1400 cm -3 to match the optical line fluxes. 


3. MODELLING THE EMISSION LINE FLUXES 

As well as aiming to fit the nebular infrared photometric 
fluxes due to dust emission, we also fitted the emission line 
fluxes from the ionized gas, principally the optical line fluxes 
measured for the entire nebula by |Smith| (12003) 1, which we 
dereddened using E(B-V) = 0.52 (Miller 1973} and the Galac¬ 
tic reddening law of Howarth ( 1983| >. We assumed an intrin¬ 
sic Case B Ha/H/3 flux ratio of 2.85 in order to determine the 
[N II] 6584,6548 A contribution to the dereddened combined 
Ha+[N II] flux, and a Case B H 7 /H /3 flux ratio of 0.47 in order 
to determine the [O III] 4363 A contribution to the dereddened 
H7+[0 III] flux. To diagnose the abundances of carbon and 
argon, lines of which did not fall within the spectral coverage 
of Smith (2003), we fitted the [C 1 ] 9824, 9850 A lines and 


the [Ar III] 7136 A line, using the dereddened line intensi¬ 


ties relative to H/3 measured by Rudy et al. (1994) for Knot 6 
(FK 6) o f|Fesen & Kirshner| ( 1982} . We note that for FK 10 


Rudy et al. (1994 1 measured [C 1 ] and [Ar III] intensities rel¬ 
ative to H/3 that were 4.0 and 2.2 times higher, respectively, 
than for FK 6. We used their FK 6 relative line intensities be¬ 
cause at shorter wavelengths the FK 6 relative line intensities 
of Henry et al. ( |1984| ) show a better match to those measured 
for the entire nebula by Smith (2003 1 


As initial nebular abundances, we used the Crab Neb¬ 
ula ‘Domain 2’ heavy element abundances from Table 2 of 
MacAlpine & Satterfield (2008 1 . Adopting a distance of 

2 kpc, we fitted the dereddened total H/3 flux by varying the 

value of the density of hydrogen in the smooth shell mod¬ 

els, or within the clumps in the clumped shell models. The 

heavy element abundances were iteratively adjusted in order 

to match the observed line fluxes, including those sensitive to 

the nebular temperature. The inferred heavy element abun¬ 
dances, by number, are listed in TablejT] Table [2]presents the 
dereddened integrated nebular line fluxes, together with the 
predicted line fluxes from one smooth model (Model III) and 
from one clumped model (Model VI). The other two smooth 
models yielded line intensity results that were very similar to 
those from Model III, while the two other clumped models 
gave line intensities that were very similar to those listed for 
Model VI (and Model III). This make it clear that the spec¬ 
tral distribution of the ionizing diffuse radiation field from the 
PWN is the most important factor in determining the emitted 
nebular spectrum. 


3.1. Results from the nebular gas-phase modelling 

From Table [3] the total nebular gas mass required to match 
the observed line fluxes ranges from 15.5 M© to 49 M© 
for the three smoothly distributed models, whereas for the 
clumped models the total gas mass is only 6.9 M© to 7.0 M©. 
The clumped model gas masses are consistent with the 8- 

10 M© mass estimated for the Crab Nebula’s progenitor star 
(Smith 2013 and references therein), whereas the smooth 
models are clearly ruled out. Optical emission line images of 
the Crab Nebula (e.g. Figure[T|i also clearly rule out a smooth 
distribution for the emitting gas. 


Table[4]presents the global elemental ion fractions obtained 
from our clumped Model VI for the Crab Nebula (the ion frac¬ 
tion patterns are very similar for its smooth model equivalent. 
Model III). Most elements are found to have a neutral fac¬ 
tion of about 10%, with the exception of helium, whose neu¬ 
tral fraction is significantly higher, at 33%. A consequence 
of the high helium neutral fraction in the Crab is that stan¬ 
dard abundance analyses based on recombination lines of H + , 
He + and He 2+ will underestimate the true He/H ratios. We 
find a helium mass fraction of 85% (Table [3), in agreement 
with the 89% derived by MacAlpine & Satterfield (2008 1 from 
their photoionization modelling of spectra from many loca¬ 
tions within the nebula. They also found that the majority of 
their locations (their ‘Domain 2’) had C/O ratios greater than 
unity, both by number and by mass. Our clumped Models V 
and VI for the entire nebula are consistent with those results, 
yielding a C/O ratio of 1.65 by number. The mass ratio of 
C/(H+He) is enhanced by a factor 6.2 in the Crab Nebula rel¬ 
ative to the solar abundances of Asplund et al. ( |2009| ), while 
the 0/(H+He) mass ratio is enhanced by only a factor of 2.3. 
The corresponding mass ratios of neon, sulphur and argon for 
the Crab are enhanced by factors of 3.8, 4.9 and 3.1 relative 
to solar, while nitrogen is depleted by a factor of 1.7. 

C/O mass ratios exceeding unity are not currently predicted 
by any supernova nucleosynthesis models, for any mass of 
progenitor star. For the CCSN yields tabulated by |Woosley| 
& Weaver ([1995]), the ejecta C/O mass ratio did increase with 


decreasing progenitor mass but for the lowest mass cases that 
they treated (11-12 M 0 ) the predicted C/O mass ratio was 
0.39, for the case of initial solar metallicity, and the predicted 
carbon yield was only 0.053 M 0 , i.e. lower than the Crab 
Nebula’s gas-phase carbon mass alone of 0.099 M 0 . In addi¬ 
tion, their 11 M 0 model predicted an ejecta He/H mass ratio 
of 0.67, versus the very much larger He/H mass ratio of 7.3 
found here for the Crab Nebula. For the lowest progenitor 
mass model (13 M 0 ) of Thielemann et al. (|1996|), an even 
lower ejecta carbon mass and C/O mass ratio was predicted 
than for the Woosley & Weaver (1995 1 model of the same 
mass. The 13 M 0 model of Nomoto et al. ( 2006| ) predicted 
a He/H mass ratio of 0.7 and a C/O mass ratio of 0.5, also 
too low compared to Crab Nebula ratios. The trend for pre¬ 
dicted carbon yields to increase with decreasing progenitor 
mass suggests that it would be useful to calculate yields for 
CCSN progenitor masses down to 8 M 0 . However, we con¬ 
clude that since no existing CCSN yield predictions match the 
case of the Crab Nebula, they therefore do not provide useful 
constraints on the total masses of heavy elements that could 
be in the gas phase or tied up within dust grains in the Crab 
Nebula. 

From an empirical analysis, Fesen et ak]( 1997j > estimated 
a total gas mass of 4.6±1.8 M 0 for the Crab Nebula, of 
which 1.5 M 0 was estimated to be in neutral filaments. From 
our clumped photoionization model we find a total gas mass 


7.0 M 0 , of which 2.1 M 0 is neutral (Tables 3 and 4). Fesen 
et al.l(jT997j( used an observed total H/3 flux of 1.78± 0.2(J / 
10~ n ergs cm -2 s _1 , from MacAlpine & Uomoto (1991), 
versus the va lue of 1.38± 0.07 x IQ -11 erg s cm~ z s 1 from 
|Smith| ( |2003] > used here. |Kirshner| ( 11974] ) measured a total 
H/3 flux of 1.30±0.40 x IQ -11 ergs cm~^ s -1 , ' ^ rrrn 


while 


David- 


(1987) estimated 1.16±0.12 x 10 11 ergs cm 2 s x . We 


adopt a factor of 1.15 uncertainty in the total H/3 flux, which 
corresponds to a factor of 1.15 1 / 2 * * * * * * * 10 = 1.07 uncertainty in the 
total nebular gas mass of 7.0±0.5 M 0 . 
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TABLE 2 

Dereddened and modelled absolute H/3 fluxes, plus dereddened and modelled line strengths relative to H/3. 


Species 

Wavelength 

[A] 

Dereddened 

Flux 1 

Modelled Flux 
Smooth III 

Dered/Model 

Modelled Flux 
Clumped VI 

Dered/Model 

H/3 

4861 

7.85x 10 -11 

6.32 x 10 -11 

1.24 

7.24 x 10 -11 

1.08 

[O ll] 

3726+3729 

18.11 

20.1 

0.90 

18.86 

1.07 

[Ne III] 

3869 

4.65 

3.79 

1.23 

3.99 

1.17 

[S II] 

4069+4076 

0.37 

0.32 

1.16 

0.36 

1.03 

[O III] 

4363 

0.57 

0.54 

1.06 

0.47 

1.20 

He I 

4471 

0.37 

0.43 

0.86 

0.37 

1.01 

He II 

4686 

0.78 

0.79 

0.98 

0.79 

0.99 

H/3 

4861 

1.00 

1.00 

1.00 

1.00 

1.00 

[O III] 

5007 

11.92 

9.57 

1.24 

9.9 

1.19 

[Nl] 

5198+5200 

0.13 

0.14 

0.93 

0.15 

0.87 

[Nil] 

5755 

0.093 

0.086 

1.08 

0.076 

1.22 

[O I]+[S III] 

6300,6363+6312 

1.23 

1.63 

0.75 

1.21 

1.02 

Ha 

6563 

2.85 

2.92 

0.98 

2.95 

0.97 

[Nil] 

6548+6584 

6.87 

6.38 

1.08 

5.70 

1.06 

[S II] 

6717+6731 

4.31 

3.98 

0.90 

4.08 

0.94 

[Ar III] 

7136 

0.34 

0.33 

1.04 

0.41 

0.84 

[Cl] 

9824+9850 

0.36 

0.66 

0.55 

0.28 

1.29 


1 Integrated line fluxes for entire nebula are fromjSmith (2003 [ , dereddened using E(B-V) = 0.52; except for [C 1] 9824+9850 and [Ar III] 7136 relative fluxes, which are from 
|Rudy et ai~|| 1 99+}. First row fluxes are in ergs cm— s — ; the fluxes in the remaining rows are relative to H/3=1.00. 


TABLE 3 

Gas phase elemental masses in the Crab Nebula 


? 

Model I 
Mass (Mg) 

Model II 
Mass (Mg) 

Model III 
Mass (Mg) 

Model IV 
Mass (Mg) 

Model V 
Mass (Mg) 

Model VI 
Mass (Mg) 

Hydrogen 

1.8 

5.53 

4.47 

0.8 

0.81 

0.81 

Helium 

13.2 

42 

33.97 

5.86 

5.99 

5.99 

Carbon 

0.21 

0.62 

0.49 

9.3xl0 -2 

9.91 xl0“ 2 

9.91xl0 -2 

Nitrogen 

6.3xl0 -3 

1.2xl0~ 2 

9.39x 10 -3 

2.8x 10 -3 

2.84x 10 -3 

2.84xl0 -3 

Oxygen 

0.2 

0.71 

0.5 

9.2x 10 -2 

7.94x 10 -2 

7.94xl0 -2 

Neon 

0.05 

5.0x 10 -2 

0.04 

3.2x 10 -2 

7.8x 10 -2 

7.8x 10 -2 

Sulphur 

0.02 

8.8x 10 -3 

7.15x 10 -3 

1.02x 10 -2 

1.04x 10 -3 

1.04xl0 -3 

Argon 

3.1xl0 -3 

1.2 —3 

6.48x 10 -3 

1.4x 10 -3 

2.9x 10 -3 

2.9x 10 -3 

Total 

15.5 

48.9 

40.1 

6.89 

7.02 

7.02 


TABLE 4 

Gas-phase elemental ion fractions for the best-fit clumpy model vi 


Neutral 1+ 2+ 3+ 4+ 5+ 


Hydrogen 

0.130 

0.870 


Helium 

0.332 

0.630 

3.77x10 

Carbon 

l.OlxlO -2 

0.730 

0.248 

Nitrogen 

1.04x 10 -2 

0.708 

0.237 

Oxygen 

0.144 

0.721 

0.107 

Neon 

0.114 

0.772 

0.113 

Sulphur 

0.198 

0.440 

0.299 

Argon 

2.31xl0 -5 

0.116 

0.702 


2.08x 10 -2 

2.01x10-® 

1.04x 10 -10 

5.39x 10 -3 

1.17X10- 6 

2.34x 10 -9 

2.75x 10 -3 

1.05x10-8 

7.37x10-® 

3.72x 10 -4 

4.10x10-8 

9.93x 10 -9 

7.05x 10 -3 

3.34xl0- 5 

5.66x10-® 

0.178 

2.31X10 -3 

4.25X10- 5 


The photoionization models described above included the 
dust grain components that are described in the next Section. 
However, running the models without dust made an insignifi¬ 
cant difference to the emission line fits, i.e. dust does not com¬ 
pete significantly for the photons that determine the global 
gas-phase ionization balance and line emission from the neb¬ 
ula. By contrast, for a model run without gas, the dust emis¬ 
sion was a factor of 1.17 lower than for the gas+dust model, 
indicating that absorption of nebular gaseous emission lines 
makes a significant contribution to the dust luminosity. 


3.1.1. Argon and ArH + 


Barlow et al. (2013) discovered the noble gas molecular ion 
36 ArH + in the Crab Nebula, via the detection of its J=1 -0 and 
2-1 rotational emission lines in Herschel -SPIRE FTS spectra. 
We therefore included argon in our photoionization modelling 
of the nebula. As noted above, for Knot FK 6, typical of the 
nebula as a whole, argon’s mass fraction was found to be en¬ 
hanced by a factor of three relative to solar. However, [Ar III] 
relative line intensities at Knot FK 10, where ArH + emission 
is strongest, are a factor of two higher than at FK 6, suggest¬ 
ing a larger enhancement of the argon abundance there. Fol- 
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lowing the detection of ArH + emission in the Crab Nebula, 
Schilke et al.|( |2014] ) were able to use a previously unidenti¬ 
fied interstellar absorption feature, now identified as due to 
ground-state absorption by the J= 1 -0 rotational line of ArH + , 
to diagnose the physical conditions in the absorbing interstel¬ 
lar clouds. They concluded that the formation reaction H 2 
+ Ar + —► H + ArH + must take place in regions where hy¬ 
drogen is overwhelmingly in the form of neutral atoms. This 
is relevant to the Crab Nebula, since photoionization models 
predict the existence of significant zones of atomic hydrogen 
(see Table 4 and [Richardson et al.| ( [2013| ). Richardson et al. 
concluded from their models that the Crab Nebula knots from 
which H 2 line emission had been detected were almost en¬ 
tirely atomic. The likelihood that ArH + also forms and emits 
in these regions should be investigated with further modelling. 

4. MODELLING THE DUST COMPONENT 

Our smooth and clumped radiative transfer models treat 
both gas and dust and, as described below, have been run with 
a wide range of dust grain parameters in order to find optimum 
fits to the observed infared spectral energy distribution of the 
Crab Nebula in order to diagnose the mass of dust present. 

From an analysis of Spitzer spectra, Temim et al. (2012 1 
found the majority of the warmer dust in the Crab nebula to 
be located in the clumpy filamentary structures. This con¬ 
clusion was supported by synchrotron-subtracted Spitzer and 
Herschel images presented by Gomez et al. ( 2012] ), which 
showed both the warm and cool dust to be concentrated in 
the nebular filaments. Figure [I] shows the far-infrared dust 
emission from the Crab Nebula to be closely aligned with 
the knots and filaments that dominate optical emission line 
images of the nebula. Given this evidence and the fact that 
smooth models require an implausibly large nebular gas mass, 
we will concentrate below on the results from our clumped 
gas+dust models. 

4.1. The grain species and their optical constants 

As discussed in Section 3.1, the Crab Nebula is carbon-rich, 
with C/0>1 by number (Table [2]), though with a few oxygen- 
rich zones ( MacAlpine & Satterfield[2008| i. Since its Spitzer 
IRS spectra sBow no features at 10 or 20 pm attributable to 
silicate Si-O stretching or bending modes ( jTemim et al. 12006} 
2012 | l, we have focused largely on amorphous carbon as the 
dominant grain species. However, we did construct some 
models that included silicates, using the silicate optical con¬ 
stants of Draine & Lee 1984k, from whom we also adopted 


the optical constants for our graphite grain models. 


For their amorphous carbon model s, |Temim & Dwek 
(20131 used optical constants from Rouleau & Martin G991|l 


(their ‘AC1’) and from |Zubko et al.| ( |1996| > (‘BE’). For com 
parison purposes we also ran models using the |Rouleau| 
& Martin ( |l99lj > AC1 amorphous carbon constants, based 


largely on optical constants measured by Bussoletti et al¬ 
i i 1987 1, as well as models with the |Zubko et al.| ( | 1996[ > BE 


amorphous carbon optical constants. The latter were based 
on data measured by Colangeli et al. ( j 1995[ > for carbon par¬ 
ticles produced from burning benzene samples. We addition¬ 
ally ran models using the |Zubko et al.| ( [ 1996k AGAR’ opti- 
cal constants, based also on measurements by|Colangeli et al. 


(19951, this time for particles produced via an electrical dis¬ 


charge through carbon electrodes in argon gas. 

The ACAR and BE amorphous carbon optical constants 
of Zubko et al. ( |1996| have no data points for wavelengths 
shorter than 40 nm and 54 nm, respectively. Since a signif¬ 


icant fraction of the Crab PWN luminosity is emitted short- 
wards of these wavelengths, we extended the BE and the 
ACAR optical constants down to shorter wavelengths using 
the 2.8 nm to 30 nm amorphous carbon optical constant mea¬ 
surements of |Uspenskii et al.] ( |2006| ), which are presented in 
the Appendix. Also listed there are n and k optical constants 
suitable for |Zubko et al.| ( jj 996| BE and ACAR grains over 
the 0.35-54 nm wavelength range, obtained as described in 
the Appendix. Figure [2] shows a comparison from 0.35 nm to 
1000 pm between the absorption efficiencies of 0.1 -pm radius 
amorphous carbon grains for the supplemented Zubko et al.| 
( |1996[ ) BE and ACAR optical constants, as well as for the op¬ 
tical constants of Rouleau & Martin] ( |1991[ ). As the latter do 
not have any data points longwards of 300 pm, we extrapo¬ 
lated them to 1000 ^m by fitting power laws to their n and k 
data points from 10-300 pm, since they change smoothly over 
this range. 

Inspection of the absorption efficiencies plotted in Fig- 
ure|2]shows s ignific ant differences between the supplemented 


Zubko et al. 


Rouleau & Martin 


((1996) BE and ACAR efficiencies and those of 
(jl99l]l, especially at wavelengths below 


20 nm and longwards of 310 nm. For wavelengths below 
100 nm we prefer the supplemented BE and ACAR data, since 
we extended these below 30 nm by using experimental opti¬ 
cal constants for amorphous carbon measured by Uspenskil] 
|et al.| ( |2006[ >. In particular, the |Uspenskii et al.| ( |2006| > data 
show a much smaller discontinuity at the carbon atom K-edge 
at 282 eV (4.4 nm) than the data of |Rouleau & Martin] ( |1 99 1| ). 
Since K-shell edges correspond to the ejection by photons of 
inner shell electrons from atoms, the vast majority of the pho¬ 
ton energy does not go into grain heating but in to raising 
the K-shell electron out of its potential well. Therefore for 
grain heating calculations, the inclusion of K-shell absorp¬ 
tion peaks will significantly overestimate the amount of grain 
heating that results. 

For wavelengths longwards of 310 nm, the ‘AC1’ amor¬ 


phous carbon optical constants presented by Rouleau & Mar- 
|tin|( 1991 1 made use of laboratory measurements of extinction 


efficiencies published by Bussoletti et al. ( 1987| >. The latter 
group subsequently obtained new laboratory measurements of 
mass extinction coefficients for different types of amorphous 
carbon particles (Colangeli et al. |1995) . They noted that their 
new data agreed with the measurements of |Koike et al.|(|1980|l 
for similar particles but not with their own (Bussoletti et al. 


et al. 


1987 earlier measurements. The newer data of (Colangeli 


1995) were used to produce the amorphous carbon op¬ 


tical constants presented by Zubk o~et al.| ( |1996| l, and overall 
we consider these, and their extensions here to shorter wave¬ 
lengths, to provide the most reliable data available for amor¬ 
phous carbon. 

For our modelling we adopted a mass density of 
1.85 g cm -3 for amorphous carbon, 2.2 g cm -3 for graphite, 
and 3.3 g cm^ 3 for silicate. 

4.2. Fitting the observed infrared and submillimeter 
photometric continuum fluxes 

The model dust SEDs were fitted to the observed in¬ 
frared and submm photometric fluxes, using observations 
made by Spitzer (Temim et al. |2006|l, Hers chel ( Gomez | 
et al. 112012} and Planck (Planck Collaboration 2011 , along 
with the mean synchrotron-subtracted Spitzer- IRS contin¬ 
uum spectrum of Temim & Dwek ( 2013] ), to better constrain 
the warm dust emission. The 24-, 70- and 100-/im points 
have been corrected for line emission, using the line con- 
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Wavelength (|im) 

Fig. 2.— Qabs(A) versus wavelength for carbon grains (apart from the Draine & Lee silicate) of radius 0.1 /zm, with optical constants taken from the labelled 
sources discussed in Section 4.1. Dashed or dotted portions use extrapolated or interpolated optical constants (see text). 


TABLE 5 

Continuum IR Fluxes from the Crab Nebula 


Wavelength 

O m ) 

Total Flux 1 

(Jy) 

Uncertainty 

(Jy) 

Dust Flux 

(Jy) 

Instrument 2 

3.6 

12.6 

0.22 


Spitzer 

4.5 

14.4 

0.26 


Spitzer 

5.8 

16.8 

0.10 


Spitzer 

8.0 

18.3 

0.13 


Spitzer 

24 

46.4 

8.0 

17.2 

Spitzer 

70 

202.4 

20 

156.8 

Herschel 

too 

196.5 

20 

143.6 

Herschel 

160 

141.8 

15 

77.5 

Herschel 

250 

103.4 

7.2 

25.9 

Herschel 

350 

102.4 

7.2 

13.2 

Herschel 

350 

99.3 

2.1 

10.1 

Planck 

500 

129.0 

9.0 


Herschel 

550 

117.7 

2.1 


Planck 

850 

128.6 

3.1 


Planck 

1382 

147.2 

3.1 


Planck 


The 24-, 70- and 100-/zm fluxes have been corrected for line emission following 


Table 2 ofjGomez et al.112012 


z Spitzer data: Temim et al. 

2006 1 ; Herschel data: |Gomez et al.|(2012); Planck 

data:|Planck Collaboration 1(2011 1 . 


tinuum fluxes between 70 and 160-/an and the uncertainties 
in fitting the SED for each grain/nebular model. 

When models were initially run with a standard MRN grain 
size distribution ( [Mathis. Rumpl, & NordsieckJ 1977| >, i.e. 
n(a) oc a~ a with a = 3.5, a min = 0.005 ftm and a max = 
0.25 /im, the dust energy distribution was found to peak at 
too short a wavelength. To better fit the observed peak, which 
is at about 70 /im, the maximum grain radius had to be in¬ 
creased to provide larger, cooler, grains, and the power-law 
slope a had to be decreased, increasing the relative number of 
larger grains. Since colder dust emits less efficiently at a given 
wavelength than warmer dust, the dust mass required to fit the 
observed fluxes also had to increase. For a number of differ¬ 
ent grain types, a = 2.7-3.0 was found to provide the best fit 
to the observed SED (see Tables [6] and |7J. For values of a < 
4, the largest grains dominate the total dust mass. There is a 
degeneracy between the maximum grain size and the slope of 
the power law, a, however better fits to both the mid and far 
infrared components of the SED were achieved with a lower 
a max and a. The value of a max was varied between 0.1 and 2 
/rm, a m i n was varied between 0.0005 and 0.1 p, m and a was 
varied between 2.4 and 4. 


tribution factors listed in Table 2 of Gomez et al. (2012). 
The dust+synchrotron continuum fluxes are listed in Table [5] 
along with 24-350-/mi dust continuum fluxes obtained by 
subtracting the synchrotron fluxes listed in Table 4 of Gomez 


et al. (2012)). For a distance of 2 kpc, the total luminosity emit¬ 


ted by dust at infrared wavelengths is 1190 Lq, corresponding 
to the absorption and reradiation of 28% of the luminosity of 
the pulsar wind nebula emitted between 0.1 nm and 1.0 /im. 

Fitting the models to the observations was done by assum¬ 
ing that the uncertainties associated with each of the observed 
fluxes were Gaussian, sampling randomly within the allowed 
observational uncertainty range to generate 1000 separate ver¬ 
sions of the SED. These were compared to the model SEDs 
and the set of model parameters generating the lowest mean 
X 2 value was taken to be the most likely. The best fit models 
to the Crab Nebula’s infrared SED are shown in Figure[3]for 
a number of different grain types, while the dust parameters 
used to obtain these fits are listed in Tables [6] and [7] The un¬ 
certainties listed for the derived dust masses are based on the 
combination in quadrature of the uncertainties in the dust con- 


4.3. The Mass of Dust 

The flux from optically thin dust emission is linearly pro¬ 
portional to the total number of dust grains, irrespective of 
whether they are in clumps. The reason that our clumped 
models have larger dust masses than our smooth models (a 
factor of ~ 1.7 larger in the case of amorphous carbon mod¬ 
els III vs. VI with Zubko ACAR and BE grain constants) is 
that the short wavelength radiation that heats the grains is at¬ 
tenuated by the gas and dust within the clumps, leading to 
cooler grains in the clump interiors than would otherwise be 
the case. Since cooler grains emit less efficiently, a larger 
mass of grains needs to be accommodated when matching a 
given far-infrared flux. 

Recombination and forbidden line emissivities are propor¬ 
tional to gas density squared, enabling the higher gas density 
clumped model to fit the observed line fluxes with a plausible 
nebular gas mass while ruling out smoothly distributed mod¬ 
els because they require an implausible 16-49 M 0 of gas in 
the Crab, versus 7 M 0 of gas for the best-fit clumped models 
(Table 3). We therefore consider that only the clumped mod- 
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TABLE 6 

Dust masses for the best fit gas+dust smooth models for the Crab Nebula 


Model 1 

Cadez et al. 

2004 

shell: 0.1 pc thick at 0.55 pc radius with central heating source 


Optical Constants 

&min 


d-max 

a 

M-dust 

x 2 

Zubko et al. ACAR 

Zubko et al. BE 

Rouleau & Martin AC 1 
Mixed Chemistry 

Draine & Lee Silicate 
Draine & Lee Graphite 

0.005to o^ °- 7 ± 0 01 F™ 

0 . 005 to.ooi Mm 0.5 ± 0.01 Atm 

0.01 zb 0.01 fim 0.8 zb 0.01 fim 

0.01 fim 1.0 fim 

0.01 d= 0.01 fim 0.9 zb 0.01 ^m 

0.001 =t 0.001 fim 0.25 zb 0.01 /xm 

2.7 ± 0.1 

2.7 ± 0.1 

2.9 ±0.1 
3.0 ±0.1 

3.5 ±0.1 

2.8 ± 0.1 

0.21 ± 0.02 Mg 

0.18 ± 0.01 Mg 

0.10 ±0.01 Mg 

0.25 ± 0.02 Mg 

0.05 ± 0.01 Mg Zubko BE 
0.14 ± 0.01 Mg DL silicates 
0.33 ± 0.04 Mg 

0.11 ± 0.01 Mg 

5.54 

3.39 

5.21 

5.23 

6.11 

4.57 


Model II 

Davidson & Fesen 119851: 

shell: 2.1xl.4pc to 4.0x2.9 pc 


Optical Constants 

&min 


&max 

Q 

M-dust 

Reduced x 2 

Zubko et al. ACAR 

0.01 ± 0.01 

fim 

1.0 zb 0.01 fim 

2.9 ±0.1 

0.18 ±0.03 Mg 

9.9 

Zubko et al BE 

0.01 ± 0.01 

fim 

0.5 zb 0.01 ^m 

2.9 ±0.1 

0.14 ±0.02 Mg 

9.7 

Rouleau & Martin AC 1 

0.01 ± 0.0 

fim 

1.0 zb 0.01 fim 

3.0 ±0.1 

0.08 ± 0.01 Mg 

12.1 

Mixed Chemistry 

0.01 ± 0.0 

fim 

0.8 zb 0.01 fim 

3.0 ±0.1 

0.29 ± 0.02 Mg 

0.10 ± 0.01 Mg Zubko BE 

0.19 ± 0.01 Mg DL silicates 

6.31 

Draine & Lee Silicate 

0.01 ± 0.01 

fim 

1.0 zb 0.01 fim 

3.5 ±0.1 

0.48 ± O.lMg 

11.3 

Draine & Lee Graphite 

0.001 ± 0.001 Atm 0.25 ± 0.01 Atm 

3.0 ±0.1 

0.09 ± 0.01 Mg 

11.0 


Optical Constants 


Model ill 


Lawrence et al. 


1995 


shell: 2.3x1.7 pc to 4.0x2.9 pc 

M dust 


Reduced % 2 


Zubko et al. ACAR 

0.01 ± 0.01 Atm 

0.7 zb 0.01 fim 

2.9 ±0.1 

0.14 ±0.04 Mg 

5.22 

Zubko et al. BE 

0.005 ± 0.005 Atm 

0.5 zb 0.01 fim 

2.8 ± 0.1 

0.11 ± 0.02 Mg 

5.97 

Rouleau & Martin AC 1 

0.01 ± 0.01 fim 

0.7 zb 0.01 fim 

3.0 ±0.1 

0.06 ± 0.01 Mg 

4.89 

Mixed Chemistry 

0.01 fim 

0.8 fim 

3.0 ±0.1 

0.21 ± 0.02 Mg 

0.07 ± 0.01 Mg Zubko BE 

0.14 ± 0.01 Mg DL silicates 

6.92 

Draine & Lee Silicate 

0.001 ± 0.001 Atm 

0.9 ± 0.01 A»m 

3.5 ±0.1 

0.37 ± 0.06 Mg 

5.14 

Draine & Lee Graphite 

0.001 ± 0.001 Atm 

0.25 ± 0.01 Atm 

2.9 ±0.1 

0.07 ± 0.01 Mg 

6.66 



Fig. 3.— The best fit overall SEDs for clumped model VI. The SEDs corresponds to the parameters in Table [7] for the different grain types described in 
Section 4.1. The observational data points are described in Section 4.2. 





















Dust and Gas Content of the Crab Nebula 


9 


TABLE 7 

Dust masses for the best fit gas+dust clumped models for the Crab Nebula 


Model IV - 

Clumps beyond the ionising radiation: 3.0x2.0 

DC to 4.0 x 2.9 pc 


Optical Constants 

d-min 


&max 

a 


x 2 

Zubko AC 

0.07 ± 0.01 Atm 

1.0 zb 0.01 fim 

2.9 ±0.1 

0.40 ± 0.08 Mq 

5.39 

Zubko BE 

0.07 ± 0.01 /rm 

0.2 zb 0.01 fim 

2.9 ±0.1 

0.30 ± 0.06 Mq 

5.52 

Rouleau & Martin AC 

0.07 ± 0.01 Atm 

1.0 zb 0.01 fim 

3.0 ±0.10 

0.24 Mq 

4.35 

Mixed Chemistry 

0.07 ± 0.01 Atm 

1.0 zb 0.01 fim 

3.0 ±0.1 

0.78 M 0 

6.61 






0.18 zb 0.03 M 0 Zubko BE 
0.60 zb 0.03 M 0 DL silicates 


Draine & Lee Silicate 

0.07 ±0.01 Atm 

1.0 zbO.Ol ^m 

3.5 ±0.1 

1.5 Mq 

4.38 

Draine & Lee Graphite 

0.001 ±0.01 Atm 

0.25 ±0.01 Atm 

3.0 ±0.1 

0.40 Mq 

3.22 

V Clumped|Lawrence et al.|(l995) shell: 2.3x1.7 

pc to 4.0x2.9 

pc - 1.lxl.l source 


Optical Constants 

d-min 


3-max 

a 

M-dust 

X 2 

Zubko et al. ACAR 

0.005 

"*" U-UU5 um 

0.7 zb 0.01 fim 

2.7 ± 0.1 

0.27 ± 0.04 Mq 

6.08 

Zubko et al. BE 

0.005t““““ /im 

0.5 zb 0.01 fim 

2.7 ±0.1 

0.20 ± 0.03 M 0 

5.99 

Rouleau & Martin AC1 

0.01 ± 0.01 Atm 

0.8 zb 0.01 fim 

2.9 ±0.1 

0.17 ±0.03 M 0 

4.98 

Mixed Chemistry 

0.01 ± 0.01 Atm 

1.0 zb 0.01 fim 

3.0 ±0.1 

0.58 ± 0.05 Mq 

5.76 






0.13 zb 0.02 Mq Zubko BE 
0.47 zb 0.03 Mq DL silicates 


Draine & Lee Silicate 

0.01 zb 0.005 fim 

0.9 zb 0.01 fim 

3.5 ±0.1 

1.10 ± 0.19 Mq 

5.44 

Draine & Lee Graphite 

0.001 ± 0.001 Atm 

0.25 ± 0.01 Atm 

2.8 ±0.1 

0.2-± 0.03 Mq 

6.03 

Model VI Clumped 

Lawrence et al. 

1995p shell: 2.3x1.7 pc to 4.0x2.9 pc - full nebula source 

Optical Constants 

d-min 


3-max 

a 

M-dust 

x 2 

Zubko et al. ACAR 

0.005 

-(-U.UU5 

-0.001 7™ 

0.7 zb 0.01 fim 

2.7 ±0.1 

0.25 ± 0.04 M 0 

5.72 

Zubko et al. BE 

0.005 

-|-0.005 

-0.001 

0.5 zb 0.01 fim 

2.7 ±0.1 

0.18 ±0.03 Mq 

4.87 

Rouleau & Martin AC1 

0.01 ± 0.01 Atm 

0.8 zb 0.01 fim 

2.9 ±0.1 

0.15 ±0.03 M 0 

4.38 

Mixed Chemistry 

0.01 ± 0.01 Atm 

1.0 zb 0.01 fim 

3.0 ±0.1 

0.50 ± 0.05 Mq 

6.62 






0.11 ± 0.02 Mq Zubko BE 
0.39 ± 0.03 Mq DL silicates 


Draine & Lee Silicate 

0.01 ± 0.01 Atm 

0.9 zb 0.01 fim 

3.5 ± 0.1 

0.98 ± 0.19 Mq 

5.12 

Draine & Lee Graphite 

0.001 ± 0.001 Atm 

0.25 ± 0.01 Atm 

2.8 ±0.1 

0.17 ±0.03 Mq 

6.42 
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Wavelength (urn) 

Fig. 4.— The best fit synchrotron subtracted clumped models V (left column) and VI (right column) dust fluxes. The dust SEDs corresponding to the parameters 
in Table [7] for the different grain types described in Section 4.1. The observational data points are described in Section 4.2 plotted with the photometric points 
and Spitzer-IRS spectrum. 


els in Tableware realistic. In addition, since the Crab Nebula 
has carbon-rich gas-phase abundances (C/O > 1; Table 1), 
the models with carbon grains are preferred over those with 
silicates. 

Models V and VI have the same distribution of dust and 
gas, but have different heating sources, with model V hav¬ 
ing a centrally located source 1.1 x 1.1 pc in diameter whilst 
model VI has the clumps embedded in a source that extends 
out to 4.0x2.9 pc. The extra heating caused by the clumps 
being embedded in the radiation source rather than outside 
it means that Model VI requires less dust to fit the observed 
SED than Model V. The spectral shape and luminosity of the 
radiation field have a far greater effect on the mass of dust 
derived for the Crab Nebula than any of the geometrical and 
density effects investigated. 

Focu s ing on the carbon grain model s, since the|Rouleau &| 
Martin] (19911 amorphous carbon and Draine & Lee ( 1984) 
graphite optical constants both include inappropriate K-shell 
absorption peaks (Figure [2j that in fact do not contribute sig¬ 
nificantly to grain heating (see the discussion in Section 4.1), 
the Zubko et ak ( |1996) l BE and ACAR amorphous carbon 
models are our preferred grain species, in clumped Models 
V or VI. These models yield a total dust mass in the Crab 
Nebula of 0.18 - 0.27 M Q (Table|7). 

Since MacAlpine & Satterfield (2008 i found a few O-rich 


Draine & Lee ([1984 1 silicate models, 


regions in the predominantly C-rich Crab Nebula, a further 

D, 


et al. ( 

1996 ) I 

& Lee 

'(M4) 

and 0.47 M 0 , 


IE amorphous carbon and 0.39 M 0 of Draine 
silicates for Model Vi’s geometry, or 0.13 M 0 


and 0.47 M 0 , respectively, for Model V’s geometry. In order 


to compare with our 
also ran models with silicate optical constants from |Laor &| 


Draine ( 1993) . The resulting silicate dust mass fits were 6% 


higher than those found with the silicate optical constants of 
Drai ne & Lee| ( |1984| l . 

Allowing for the 0.099 M 0 of gas-phase carbon in the neb¬ 
ula for clumped Models V or VI (Table 3), the minimum to¬ 
tal mass of carbon in the Crab Nebula is 0.28 M 0 . As dis¬ 
cussed in Section 3.1, we do not consider that the low carbon 
yields predicted by the 11-13- M 0 core-collapse SN models 
of Woosley & Weaver ( 1995[ >, Thielemann et al. (1996 > and 
Nomoto et al. (|2006|> provide a useful constraint on the mass 


of carbon that can be in dust, since their predicted C/O and 
He/H mass ratios are much smaller than those found in the 
Crab Nebula. 

For clumped Model V with Zubko et al.| ( |1996| > BE 
amorphous carbon grains (Table [7}, the gas and dust 
masses in each clump were respectively 6.08 xl0~ 3 M 0 and 
1.68 x 10 4 M 0 , for a gas to dust mass ratio of 36, and the 
V-band dust optical depth from the edge to the centre of each 
clump was Tv = 1.12. Following the detection by |Woltjer &| 
Veron-Cetty ( 1987)> of absorption attributable to dust at the 


position of a bright [O III] filament in the Crab Nebula, [Fesen| 


& Blair ( ] 1990| ) measured angular diameters ranging from 0.9 
to 4.8 arcsec for 24 ‘dark spots’ in the Crab Nebula. For com¬ 
parison, the 0.037-pc radius clumps adopted for our clumped 
models would have an angular radius of 3.8 arcsec at a dis¬ 
tance of 2 kpc. 

4.3.1. Comparison with previous dust mass estimates 
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Since cool dust emits less efficiently than warm dust, larger 
dust masses are needed to fit far-infrared fluxes than are re¬ 
quired to fit similar fluxes at shorter infrared wavelengths. So 
observations extending out to far-infrared and submillimeter 
wavelengths are often necessary in order to fully characterize 
nebular dust masses. In our comparisons below, we will focus 
on dust mass estimates made assuming carbon grains. 

Prior to the launch of Herschel , the longest infrared wave¬ 
lengths that the Crab Nebula had been observed to were the 
IRAS 12-100 pm observations of |Marsden et al.|| ( }l984j > and 
the ISO 60-170-pm plus SCUBA 850-pm observations of 
|Green et~aL| ( |2004) l. The latter’s 60- and 100-pm fluxes were 
lower by factors of 1.5 and 1.7 respectively than the fluxes 
measured with IRAS by |Marsden et al. ( 1984| ), whose 100- 
pm fl ux was within 10% of the value measure d with Her- 
schel ( |Gomez et al. 2012| >, although |Marsden et al.| ( |1984 1 
adopted and subtracted a much larger far-infrared synchrotron 
flux component than Gomez et al. ( |2012 1 , who had accu¬ 
rate Planck submillimeter and Spitzer near-mid-infrared flux 
measurements available to define the underlying synchrotron 
spectrum. Marsden et al. (1984 1 estimated a dust mass of 
10 :i M 0 for 80 K grains having a A -1 emissivity, or 0.3 M-. 
for 50 K grains having a A -2 emissivity. Temim et al. 


( |2012| ) fitted Spitzer data that extended out to 70 pm with 
(Fi)x 10 -3 Mg of 60±7 K Zubko amorphous carbon dust. 

The advent of 70-500-pm Herschel data enabled Gomez 
et al.| p012] > to fit two modified blackbodies to the Spitzer and 
Herschel infrared and submillimeter SED of the Crab Nebula. 
For the amorphous carbon case the blackbodies were modi¬ 
fied by the wavelength dependence of the absorption coeffi¬ 
cients of Zubko et al.| ( |l996] ) BE amorphous carbon, with the 
warmer (63±4 K) and cooler (34±2 K) components requir¬ 
ing 0.006±0.02 and 0.11±0.01 Mg, respectively, of carbon 
grains, i.e. the same as the 0.11 ±0.02 M 0 of BE amorphous 
carb on dust required by o ur smooth model III (Table |6|» . 

|Temim & Dwek|(j2013|) obtained a lower carbon dust mass 
by fitting the Gomez et al. (20121 infrared SED with amor¬ 
phous carbon grains having a power-law distribution of grain 
radii whose radiative equilibrium was calculated assuming 
heating by a central point source whose spectrum matched 
that of the pulsar wind nebula (PWN). For their best-fit mod¬ 
els for the SED, the grains were at a distance of 0.5-0.7 pc 
from the center of the nebula, corresponding to ~0.20-0.25 
of the nebular radius. They derived a mass of 0.02 M 0 
for |Rouleau & Martin ( 1991|> amorphous carbon grains, or 
0.04 M 0 for|Zubko et al. ([1996 ) BE amorphous carbon grains. 

Our smooth Model I aimed to mimic the geometry adopted 
by Temim & Dwek| ( 12013} but has a diffuse PWN ionizing 
radiation source instead of a centrally located point radiation 
source. We obtained dust masses of 0.10 M 0 for Rouleau 


& Martin ( 199T]l_amorphous carbon grains and 0.18 M 0 for 


Zubko et al.[ ( |1996j > BE amorphous carbon grains. The dif¬ 
ference between these two dust masses may be attributable to 
the inclusion in the |Rouleau & Maitin| ( jl99r| data of an over¬ 
large absorption cross-section for grain heating at the K-shell 
edge of atomic carbon (see Section 3.1), together with the use 
by Zubko et al. ( |1996| ) of improved optical and longer wave¬ 
length data from the Lecce group, compared to the earlier data 
from the same group used by Rouleau & Martmjl 1991 1. When 
we ran a smooth Model I with Rouleau & Martin (1991) amor¬ 
phous carbon dust whose large K-shell absorption peak (Fig¬ 
ure 2) had been replaced by an interpolation of the underlying 
absorption efficiency, the mass of dust required to fit the in¬ 


frared SED increased from 0.10 M 0 (Table [6]) to 0.13 M 0 . 

The filaments and clumps of the Crab Nebula with which 
the dust is associated extend all the way to the outer edges 
of the nebula (Figure 1), inconsistent with the 0.55-0.65-pc 
shell geometry of Model I, which also required an implausibly 
large gas mass (15.5 M 0 ; Table 3). 

5. CONCLUSIONS 

We have constructed a series of radiative transfer models 
to determine the mass of dust present in the Crab Nebula su¬ 
pernova remnant. In the preferred models the gas and dust are 
located in clumps within an ellipsoidal diffuse synchrotron ra¬ 
diation source, powered by the pulsar wind nebula. The mod¬ 
els are insensitive to the inner axis diameters from which the 
clump distributions extend. 

Models with a smooth distribution of material require 0.11- 
0.21 M 0 of|Zubko et al .|(fl 996) BE or ACAR amorphous car¬ 
bon, respectively, or 0.33-0.48 M 0 of|Draine & Lee|(1984) 
silicates, to fit the infrared and submilhmeter SED defined 
by the Herschel and Spitzer observations of the nebula. This 
compares with the 0.12±0.02 M 0 of Zubko BE amorphous 
carbon, or the 0.24^'^ M 0 of Weingartner & Draine (2001) 


silicate dust, derived by Gomez et al. ( |2012[ > from two- 
component blackbody fits modified by the mass absorption 
coefficents for those materials. 

Our smooth distribution models required implausibly large 
nebular gas masses of 16-49 M 0 to fit the integrated optical 
line fluxes measured by Smith ( 2003| ) for the Crab Nebula, 
much larger than the 8-10 M 0 initial mass usually estimated 
for the progenitor star, whereas our clumped models for the 
gas and dust, more consistent with the filamentary appear¬ 
ance appearance of the nebula, required only 7.0±0.5 M 0 
of gas to match the integrated nebular emission line fluxes. 
The clumped model V and VI infrared SED fits, which are 
therefore preferred over those from the smooth models, re¬ 
quired either 0.18-0.20 M 0 (BE) or 0.25-0.27 M 0 (ACAR) 
of Zubko amorphous carbon, 0.98-1.10 M 0 of Draine & Lee 
silicate, or, for mixed chemistry dust, 0.11-0.13 M 0 of Zubko 
BE amorphous carbon plus 0.38-0.47 M 0 of silicates. Since 
our photoionization modelling yielded an overall gas-phase 
C/O ratio of 1.65 by number for the Crab Nebula, the clumped 
model dust masses obtained using just amorphous carbon, or 
amorphous carbon plus silicates, are favoured over silicate- 
only models. The total nebular mass (gas plus dust) is es¬ 
timated to be 7.2±0.5 M 0 . The Crab Nebula’s gas to dust 
mass ratio of 26-39 (depending on the exact grain type) is 
about 5-7 times lower than for the general ISM. As discussed 
in the Introduction, CCSN ejecta dust masses of 0.1 M 0 or 
more, a constraint satisfied by the Crab Nebula, Cas A and 
SN 1987A, can potentially make a significant contribution to 
the dust budgets of galaxies. 

Our best fit power-law grain size distributions, n(a) oc a ~ a , 
had a ~ 3, so that the majority of the dust mass resides in the 
largest particles, with a max = 0.5-1.0 pm. Larger particles 
better withstand destruction by shock sputtering, for which 
the rate of reduction of grain radius, da/dt, is independent 
of the grain radius a, so that the smallest particles disappear 
first. The preponderance of larger particles in the Crab Neb¬ 
ula’s dust, and the fact that they are in clumps, can help their 
longer-term survival when they eventually encounter the in¬ 
terstellar medium (Nozawa et a l.||2007| l. 

A mass of 8-13 M 0 has previously been estimated for the 
Crab Nebula progenitor star (Heste r[2008 Smith|201 3[ ). The 
fact that earlier nebular mass estimates have fallen well short 
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of this mass range had been used as one of the arguments that 
faster moving material must exist beyond the main nebular 
boundaries (see e.g. Hester 2008). Arguments against that 
conclusion have however been presented by Smith ( |2013| >. 
The total nebular mass of (7.2±0.5) Mq derived here, com¬ 
bined with a pulsar mass of at least 1.4 M 0 , implies a total 
mass of at least 8.6 Mq, removing a nebular mass deficit as 
an argument for the existence of additional material beyond 
the visible boundaries of the Crab Nebula. 


We thank Dr Tea Temim for comments that helped improve 
the paper and for making available the_mean synchrotron- 
subtracted Spitzer-IRS spectrum of |Temim & Dwek ( 2013[ ). 
We thank Antonia Bevan, Barbara Ercolano, Haley Gomez, 
Oskar Karczewski, Mikako Matsuura, Bruce Swinyard and 
Roger Wesson for discussions about MOCASSIN, dust and su¬ 
pernova remnants. 
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APPENDIX 

Appendix: Amorphous Carbon EUV and X-ray Optical Constants 


Table 8 lists the values of n and k measured by [Uspenskii et al.](|2006| ) between 2.8 nm and 30 nm for an amorphous carbon 
sample. It also lists extrapolated n and k values for the |Zubko et al. I fl9961 ACAR amd BE amorphous carbon samples, obtained 
by fitting power-laws to the short wavelength ends of their n and k distributions and then extrapolating these from their shortest 
wavelength points, at 40 nm and 54 nm, respectively, to shorter wavelengths until they intersected the n and k data of Uspenskii 


et al. 


et al. 


(2006), which were then used from the intersection wavelength down to 2.8 nm. Power-law extrapolations of the 


1 2006) n and k data were used from 2.8 nm down to 0.35 nm. 


Uspenskii 


Wavelength 

n 

k 

n 

k 

n 

k 

(nm) 

Uspenskii 

Uspenskii 

Zubko ACAR 

Zubko ACAR 

Zubko BE 

Zubko BE 

0.3 



9.970E-01 

1.72E-06 

9.970E-01 

1.72E-06 

0.4 



9.970E-01 

3.26E-06 

9.970E-01 

3.26E-06 

0.5 



9.970E-01 

5.36E-06 

9.970E-01 

5.36E-06 

0.6 



9.970E-01 

8.04E-06 

9.970E-01 

8.04E-06 

0.7 



9.970E-01 

1.13E-05 

9.970E-01 

1.13E-05 

0.8 



9.970E-01 

1.52E-05 

9.970E-01 

1.52E-05 

0.9 



9.970E-01 

1.98E-05 

9.970E-01 

1.98E-05 

1.5 



9.970E-01 

6.16E-05 

9.970E-01 

6.16E-05 

3.55 

9.970E-01 

4.477E-04 

9.970E-01 

4.477E-04 

9.970E-01 

4.477E-04 

3.76 

9.970E-01 

5.226E-04 

9.970E-01 

5.226E-04 

9.970E-01 

5.226E-04 

3.98 

9.971E-01 

6.004E-04 

9.971E-01 

6.004E-04 

9.971E-01 

6.004E-04 

4.13 

9.980E-01 

2.000E-03 

9.980E-01 

2.000E-03 

9.980E-01 

2.000E-03 

4.21 

9.980E-01 

1.700E-03 

9.971E-01 

1.700E-03 

9.971E-01 

1.700E-03 

4.27 

9.980E-01 

1.000E-03 

9.980E+00 

1.000E-03 

9.980E+00 

1.000E-03 

4.35 

1.000E+00 

3.000E-03 

1.000E+00 

3.000E-03 

1.000E+00 

3.000E-03 

4.42 

9.980E-01 

1.000E-04 

9.980E-01 

1.000E-04 

9.980E-01 

1.000E-04 

4.59 

9.970E-01 

1.000E-04 

9.970E-01 

1.000E-04 

9.970E-01 

1.000E-04 

4.72 

9.971E-01 

1.062E-03 

9.971E-01 

1.062E-03 

9.971E-01 

1.062E-03 

5.00 

9.970E-01 

1.062E-03 

9.970E-01 

1.062E-03 

9.970E-01 

1.062E-03 

5.29 

9.969E-01 

1.062E-03 

9.969E-01 

1.062E-03 

9.969E-01 

1.062E-03 

5.60 

9.967E-01 

1.174E-03 

9.967E-01 

1.174E-03 

9.967E-01 

1.174E-03 

5.93 

9.964E-01 

1.297E-03 

9.964E-01 

1.297E-03 

9.964E-01 

1.297E-03 

6.27 

9.960E-01 

1.433E-03 

9.960E-01 

1.433E-03 

9.960E-01 

1.433E-03 

6.64 

9.956E-01 

1.586E-03 

9.956E-01 

1.586E-03 

9.956E-01 

1.586E-03 

7.03 

9.950E-01 

1.755E-03 

9.950E-01 

1.755E-03 

9.950E-01 

1.755E-03 

7.44 

9.943E-01 

1.949E-03 

9.943E-01 

1.949E-03 

9.943E-01 

1.949E-03 

7.88 

9.935E-01 

2.174E-03 

9.935E-01 

2.174E-03 

9.935E-01 

2.174E-03 

8.34 

9.925E-01 

2.431E-03 

9.925E-01 

2.431E-03 

9.925E-01 

2.431E-03 

8.82 

9.913E-01 

2.731E-03 

9.913E-01 

2.731E-03 

9.913E-01 

2.731E-03 

9.34 

9.899E-01 

3.087E-03 

9.899E-01 

3.087E-03 

9.899E-01 

3.087E-03 

9.89 

9.882E-01 

3.511E-03 

9.882E-01 

3.511E-03 

9.882E-01 

3.51 IE-03 

10.46 

9.863E-01 

4.009E-03 

9.863E-01 

4.009E-03 

9.863E-01 

4.009E-03 

11.08 

9.841E-01 

4.612E-03 

9.841E-01 

4.612E-03 

9.841E-01 

4.612E-03 

11.73 

9.814E-01 

5.333E-03 

9.814E-01 

5.333E-03 

9.814E-01 

5.333E-03 

12.42 

9.784E-01 

6.209E-03 

9.784E-01 

6.209E-03 

9.784E-01 

6.209E-03 

13.14 

9.749E-01 

7.251E-03 

9.749E-01 

7.251E-03 

9.749E-01 

7.251E-03 

13.91 

9.709E-01 

8.519E-03 

9.709E-01 

8.519E-03 

9.709E-01 

8.519E-03 

14.72 

9.663E-01 

1.005E-02 

9.663E-01 

1.005E-02 

9.663E-01 

1.005E-02 

15.58 

9.610E-01 

1.192E-02 

9.610E-01 

1.192E-02 

9.610E-01 

1.192E-02 

16.50 

9.549E-01 

1.417E-02 

9.549E-01 

1.417E-02 

9.549E-01 

1.417E-02 

17.46 

9.480E-01 

1.691E-02 

9.480E-01 

1.691E-02 

9.480E-01 

1.691E-02 

18.49 

9.400E-01 

2.022E-02 

9.400E-01 

2.022E-02 

9.400E-01 

2.022E-02 

19.56 

9.310E-01 

2.422E-02 

9.310E-01 

2.422E-02 

9.310E-01 

2.422E-02 

20.71 

9.207E-01 

2.909E-02 

9.207E-01 

2.909E-02 

9.207E-01 

2.909E-02 

21.92 

9.091E-01 

3.496E-02 

9.091E-01 

3.496E-02 

9.091E-01 

3.496E-02 

23.21 

8.958E-01 

4.206E-02 

8.958E-01 

4.206E-02 

8.958E-01 

4.206E-02 

24.56 

8.807E-01 

5.064E-02 

8.807E-01 

5.064E-02 

8.807E-01 

5.064E-02 

26.00 

8.637E-01 

6.101E-02 

8.637E-01 

6.101E-02 

8.637E-01 

6.101E-02 

27.52 

8.443E-01 

7.350E-02 

8.443E-01 

7.350E-02 

8.443E-01 

7.350E-02 

29.13 

8.224E-01 

8.860E-02 

8.224E-01 

8.860E-02 

8.224E-01 

8.860E-02 

40.00 



9.090E-01 

7.92E-02 

8.990E-01 

9.01 IE-2 

50.00 



8.63800E-01 

1.93800E-01 

9.410E-01 

9.780E-02 

54.00 



8.60100E-01 

2.35100E-01 

9 .18300E-01 

1.26400E-01 


TABLE 8 

Values of n and k measured by Uspen skii et al~]U006[i b etween 2.8 nm and 30 nm for an amorphous carbon sample and 
EXTRAPOLATED N AND K FOR ZUBKO_ET_AL. (| 1 996|> ACAR AND BE AMORPHOUS CARBON SAMPLES 




















